Preprint 
April 2, 2013 



Fluctuation dynamo driven by sliear-bursts in convectively-driven 

magnetoliydrodynamic turbulence 

J. Pratt\ A. Busse^, and W.-C. Miiller-' 



Max-Planck-Institut fiir Plasmaphysik, 85748 Garching, Germany and Max-Planck-Institut fiir Sonnensystemforschung, 37191 

Katlenburg-Lindau, Germany 

Faculty of Engineering and the Environment, University of Southampton, UK 

Center for Astronomy and Astrophysics, TU Berlin, Hardenbergstr. 36 10623 Berlin, Germany 



April 2, 2013 



m 

o 

(N 



ABSTRACT 



cn 

p. 
I 

o 

-(— • 



Intermittent large-scale high-shear flows are found to occur frequently and spontaneously in direct numerical simulations of statisti- 
cally stationary turbulent Boussinesq magnetohydrodynamic (MHD) convection. The energetic steady-state of the system is sustained 
by convective driving of the velocity field and small-scale dynamo action. The intermittent emergence of flow structures with strong 
velocity and magnetic shearing generates magnetic energy at an elevated rate over time-scales longer than the characteristic time 
of the large-scale convective motion. The resilience of magnetic energy amplification suggests that intermittent shear-bursts are a 
significant driver of dynamo action in turbulent magnetoconvection. 
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X-ray observations reveal that turbulent convection agitates 
the outer convection layer of stars ( jGiidel et al.| 1997[ Reiners| 
& Basri| [2007 , B6hm-Vitense1 |2008) [Siinon et al'pOOst . Mea- 



surements also show that planetary magnetic fields can change 
in magnitude and orientation (,McFadden & Merrill 1995[ |Chris- 
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tensen et al.|2009[ |Stevenson|2010[|Olson et al.|20i l| i7Dynamo 
action driven by turbulent convection is accepted as the origin 
,of solar and planetary magnetic fields. Key physical processes 
involved in turbulent convection and implicated in the amplifi- 
cation of magnetic fields remain to be identified and practically 
understood ( [Zeldovich et al.|1983[|Biskamp|2003] l. Helicity and 
buoyancy remain intensely interesting to the dynamo problem 
(|Chne et al.|[2003| [CaTtaiieo et al.|p03| [Weiss & Thompson 
[2009] i: 

Because of the inherent nonlinearity of turbulent plasma 
flows, theoretical explanation of dynamo action is often ap- 
proached by mean-field theory. Comparison with three- 
dimensional numerical simulations verifies and inspires theo- 
retical models (|Moll et al.|[20n] | Wilkin et al.|[2 007; Harder & 
|Hansen|[2005| [Stanley & Gla tzmaier 2010; Tobias et_aLJ[2011|. 
This work reports on a resilient and newly-identified feature of 
characteristic dynamo action in three-dimensional convectively- 
driven homogeneously turbulent Boussinesq magnetoconvec- 
tion based on direct numerical pseudospectral simulations using 
the magnetohydrodynamic (MHD) fluid approximation ( [Chan-[ 
[drasekhar,1961j . 

Formulation of optimal boundary conditions for simulations 
of turbulent flows is delicate because boundaries strongly influ- 
ence the structure and dynamics of the flow. The commonly- 
used Rayleigh Benard boundary conditions impose a tempera- 
ture gradient to drive turbulent convection by fixing the tem- 
perature on impermeable top and bottom boundaries. For the 
Reynolds and Rayleigh numbers achievable with current numer- 
ical capabilities, the convection-cell pattern imposed on the flow 



by Rayleigh Benard boundary conditions constrains the devel- 
opment of buoyantly-driven turbulence. 

The simulation volume considered in this work is confined 
by quasi-periodic rather than Rayleigh Benard boundary condi- 
tions. The aim is to combine the conceptual simplicity of statisti- 
cal homogeneity with a physically natural convective driving of 
the turbulent flow. The only deviation from full periodicity in our 
boundary conditions is the suppression of mean flows parallel to 
gravity; this modification suppresses the macroscopic elevator 
instability discussed by Calzavari ni et al.[ ( |2006[ ). In the flow al- 
lowed by these boundary conditions in our simulation we iden- 
tify a robust universal process, the shear-burst, that efficiently 
amplifies magnetic energy at all spatial scales in convective tur- 
bulence. The simulation model we employ is idealized, but can 
be viewed as a volume in an astrophysical or geophysical con- 
vective turbulent flow that is small in comparison to the pressure 
scale height. 

Our system allows the study of a turbulent fluctuation dy- 
namo in detail since the applied boundary conditions permit 
shear-bursts on large spatial and temporal scales without enforc- 
ing a large-scale structuring of the turbulent flow. Hundreds 
of convective time-scales prove necessary for the study of the 
shear-bursts that arise spontaneously in simulations of steady- 
state convective MHD turbulence. Shear-bursts are intermittent 
and spatially localized around high-shear flows; they are driven 
primarily at multiple large length scales, which are not neces- 
sarily continuous, and vaiy between bursts. A single burst is 
not sufficient to maintain elevated magnetic energy; however, 
shear-bursts recur frequently, providing an elevated growth of 
magnetic energy over significant periods of time. We address 
the properties of shear-bursts and their importance for the under- 
standing of the fluctuation-dynamo mechanism based on obser- 
vations from high-resolution direct numerical simulations that 
span extended periods of time. 

Article number, page 1 of pi 



scaled ViB, 

scaled BjViVj 

scaled Eb 




■N^ A • / A A A A 

T I I I I I I 

270 275 280 285 290 295 300 
t/th 



00 

d 




Fig. 1. Magnetic energy, E^, is significantly amplified over long 
times due to repeated shear-bursts in (above) simulation gl and (below) 
simulation gl 1 (parameters described in TablefTl. The signature of the 
shear-burst, identified here with arrows, is clearly visible in the simul- 
taneous peaks of the magnitude of the magnetic shear tensor V,6^, and 
magnetic stretching tensor S,V,j(^, which have been normalized to their 
initial-time values for easy comparison. The magnitude of the velocity 
shear tensor VjVj closely follows the magnetic shear. The magnitude 
of these tensors is calculated as a sum of the squares of the elements. 
Shear, and therefore magnetic stretching, drive the intermittent growth 
of magnetic energy. 



The non-dimensional Boussinesq equations for MHD con- 
vection are: 



do) , 

— - V x{vxoj+ j X B) - vV cl> -V9xgr. 

ot 

dB , 

^- - VxivxB)^ nW^B 
ot 

de , 

— -1- (u ■ V)0 = kV^e - (d ■ V)Ta 
ot 



(1) 

(2) 

(3) 
(4) 



The equations include the solenoidal velocity field v, vorticity 
<x> -V XV, magnetic field B, current j - V x B. The quantity 6 
denotes the temperature fluctuation about a linear mean temper- 
ature profile To(z) where z is the direction of gravity. The mean 
temperature, To, in equation [3] provides the convective drive of 
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the system. The temperature fluctuation, 9, term in eq. [T] is 
the buoyancy force. The vector g^ is a unit vector in the direc- 
tion of gravity. These equations are solved using a pseudospec- 
tral method with an adaptive-time-integration accomplished by a 
low-storage third-order Runge Kutta scheme (Williamson 1980 1. 
Turbulent convective motion defines the characteristic time 
and length scales of the system: the large-scale buoyancy time, 
fb - (o'g'IVrol)"'^^ and temperature gradient length scale L = 
Ttf^To where r» is defined by the root-mean-square of tempera- 
ture fluctuations 6. a is the volume thermal expansion coefficient 



at constant pressure and q the gravitatio nal acceleration ( Gibert 



et al. 



2006 



f[ 



Skandera & Miiller 2009). The magnetic field is 
given in Alfvenic units, with Alfven Mach number uq/ua = 1, 
vq - L/t\,. Three dimensionless parameters appear in the equa- 
tions: V, ?), and k. They derive from the kinematic viscosity v, 
the magnetic diff'usivity t], and thermal diffusivity k and can for- 
mally be identified as the reciprocal value of classical Reynolds 
number, magnetic Reynolds number, and Peclet numbers, re- 
spectively (see TablefTlfor definitions). 

To investigate the influence of diffusivities on the shear-burst 
phenomena, parameters v, fj, and k are varied; consequently the 
simulations probe values of the Prandtl Pr = v/k and the mag- 
netic Prandtl number Prm = v/f) spread between 0.5 and 2. The 
magnetic Prandtl number has been shown to exhibit a signifi- 



cant effect on small-scale dynamo action ( Boldyrev & Cattaneo 
[2004j Schekochihin et al._2005, ); the dependence of the dynamo 
mechanism on Prandtl number has been the subject of several 
wide-ranging investigations (' Schmalzl et al.|2()02 Maron et al. 
[2004; Simitev & Busse 2005). Realistically smaU Prandtl num- 
bers cannot be reached with contemporary computer capabili- 
ties; in the solar convection zone generally expected values are 
Prm ~ 10"' - 10"^ and in the earth's core Prm ~ lO"*". Simul- 
taneously, Reynolds numbers are generally expected to be larger 
than can be computationally reached: Re ~ 10'^ in the solar 
convection zone and Re ~ 10^ in the earth's core. Because of 
this discrepancy, the dynamical ranges of fluctuations in modern 
simulations are smaller than those expected in real systems. Our 
simulations thus present a first impression of the role of shear- 
bursts for astrophysical dynamos. 

The numerical turbulence data in this work results from a set 
of simulations conducted with grid size 512-', which constitutes 
high-resolution for the extremely long times treated here. These 
simulations are performed in a quasi-periodic slab or cube; the 
cube has a side of 27r, and the slab has slightly larger x- and y- 
directions of (27r)^^^ to allow for the well-defined initial onset of 
the convective instability driving the turbulence (| Chandrasekhar| 
|1961) l. Our boundary conditions inhibit the formation of vis- 
cous boundary layers, which appear when impermeable bound- 
ary conditions are employed. The dissipative coefficients v, fj, 
and k parametrize the extent of the turbulent inertial range of 
scales, and in each simulation are chosen to be as small as possi- 
ble so that the resolution constraint of kmax^l > 1 -5 is still satisfied 
(jPope 2000). Here, fcn^ax is the largest resolvable wavenumber al- 
lowed by numerical resolution and 77 is the Kolmogorov length 
scale. The simulation parameters are summarized in Table [T] 
The Rayleigh number, characterizing the dynamical importance 
of buoyancy in Rayleigh-Benard configurations, is of limited in- 
formative value for the present quasi-periodic system. 

The initial state of the simulations contains fluctuations in 
a number of small k modes for the velocity, magnetic field, 
and temperature. The initial ratio of kinetic to magnetic en- 
ergy of turbulent fluctuations is 10^ with the kinetic energy 
of order unity. Fig. [2] shows a typical example of the ini- 
tial time-evolution taken from simulation g2 of kinetic, Ey = 
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Table 1. Dimensionless simulation parameters: magnetic Reynolds number Rem, Peclet number Pe, Prandtl number Pr = vik, magnetic Prandtl 
number Ptm = v/Zy, Rayleigh number Ra. The Reynolds number is defined as Re = (Ev + EB)^/y(ev + fe) where e represents energy dissipation. 
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Fig. 2. Total kinetic, magnetic and thermal energies during the initial 
growth phase of the dynamo in simulation g2. 



V ' J dV\/^/2, magnetic, E^ = V ^ J dVB^, and temperature, 

Ej — V ' ^ AV0^. The thermal energy should be understood 
as the variance of temperature fluctuations. Magnetic energy 
rises quickly due to small-scale dynamo action and saturates at 
EsjIEb ~ 0(1) characteristic of the quasi-stationary turbulent 
state of the MHD flow. 

In Fig. [2] the steady-state global energies of the system 
evolve with fluctuations due to the convective motion. After a 
simulation reaches steady-state, energies fluctuate on the order 
of 10%, with a period of a 1-2 buoyancy times (see also Figure 
5 of Cattaneo 1999a). During steady-state turbulent convection, 
we begin to observe a pattern of spontaneous longer periods (5- 
20 fb) of significant growth in the global magnetic, thermal and 
kinetic energies. The energy variation during one such period 
can reach 10 times the steady-state energy level. For exam- 
ple, these periods of unusually elevated energy-growth occur 15 
times over a time-span of 225 fb in the simulation gl. The sys- 
tem can be regarded as statistically steady over periods of time 
significantly longer than the duration of a shear-burst. 

A shear-burst centers around a period of increased growfli 
magnetic energy that is accompanied by growth of both magnetic 
shear and magnetic stretching. Preceding the growth of mag- 




Flg. 3. A hot (black) flow moves upward vigorously, shearing against 
a cold (light blue) flow that moves downward. The shear along these 
opposing flows drives energy production during shear-bursts. 



netic energy, a coherent flow structure forms that has the appear- 
ance of high-velocity, thermally-polarized, coherent streams, in 
contrast to the typical situation with many smoothly convect- 
ing plumes of hot and cold fluid. These high-energy streams 
are strongly aligned in space, producing regions of high and in- 
creasing velocity and magnetic shear. The nonlinear shape and 
orientation of the high-energy streams differ for each shear-burst, 
displaying no preferred direction. 

Shear causes magnetic field-line stretching, and thus ampli- 
fication of magnetic energy (Childress & Gilbert 1995 Cattaneo 
|1999b| l. In Fig. [Tjeach shear-burst can be defined by a peak in 
magnetic shear that correlates with an increase in magnetic en- 
ergy. Fig. |4] allows for closer inspection of the increase of mag- 
netic, kinetic, and thermal energies for a typical isolated shear- 
burst in simulation g8; between t - 225 and t - 235 the en- 
ergies increase by a factor of three. Individual shear-bursts can 
last from a couple buoyancy times to a couple tens of buoyancy 
times. At the peak of magnetic and kinetic energies, high energy 
hot and cold shearing becomes violent enough to break up the 
flow structures. The peak of global energy in Fig. |4]represents 
such an explosion of flow structures. The break-up of the fast 
streams spurs a slow decline in global energies. After the shear- 
burst has dissipated, the energies dissipate until the steady-state 
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Fig. 4. The global kinetic energy E^, magnetic energy Eg , and thermal 
energy Ej scaled to run between and 1, for the span of a typical shear- 
burst in simulation g8. 



level maintained by the fluctuation dynamo has been reached. 
Shear-bursts can overlap in time, and also can occur closely in 
sequence, as shown in Fig. [1] 

The lifetime and magnitude of the energy growth, in particu- 
lar the peak in thermal energy, can vary greatly between simula- 
tions and even between shear-bursts in the same simulation. This 
shows no apparent dependence on the Prandtl numbers. That 
the Prandtl numbers do not directly impact the shear-burst phe- 
nomena is surprising because the Prandtl numbers express the 
ratio of turbulent intensities and dynamic ranges of the respec- 
tive turbulent fields. This relationship between Prandtl numbers 
and the turbulence can be understood by relating the Prandtl 
numbers to the ratios of Reynolds numbers, Pr = Pe/Re and 
Prm = Rem /Re, where the Peclet number can be regarded as the 
same structure as a Reynolds number for thermal fluctuations. 

The characteristic length-scale of Boussinesq convection is 
the Bolgiano-Obukhov length, ^bo = fv /^x ^^'^^ separates 
convectively-driven scales of the flow { > ^bo from the range 
of scales where the temperature fluctuations behave as a passive 
scalar £ < ^bo- In this definition e^ and er are the kinetic and 
thermal energy dissipations respectively. Typically in our con- 
vection simulations /'bo is comparable to the system size so only 
the largest scales in the flow are convectively driven. The shear- 
bursts are large scale phenomena by this classification, but not 
dominated by convective motion. This suggests an explanation 
for the apparent insensitivity to changes in the Prandtl numbers. 

Although insensitive to the Prandtl numbers, the shear-bursts 
interact nonlinearly with the turbulent environment, mainly via 
large-scale magnetic structure. This is reflected in the behav- 
ior of magnetic helicity, Hm - V"' J dVA ■ B, which measures 



the linkage and knottedness of the magnetic field-lines (Biskamp 
|2000||Moff'att|1978| . A signature of the shear-burst is the growth 
of global magnetic helicity as the shear flows strengthen. Mag- 
netic helicity is not conserved in the dissipative system we study, 
and this growth of magnetic helicity typically exceeds more than 

Article number, page 4 ofpl 



t/tb 



Fig. 5. The characteristic double-peak from two closely spaced 

shear-bursts in simulation gll is visible in magnetic helicity. The mag- 
netic helicity Hm at low-wavenumbers k is shown scaled to its initial 
time-value with sign preserved. Magnetic helicity quantities are shown 
scaled to their initial time values in a way that preserves sign. The mag- 
nitude of the magnetic stretching tensor is shown shifted and scaled to 
fit the scale of the magnetic helicity, to provide a time-reference for the 
shear-bursts. 



a standard deviation from the average magnetic helicity over the 
time-span of the simulation. A peak of global magnetic helic- 
ity frequently shortly precedes or coincides with a shear-burst. 
Fig. [5] shows the typical time-evolution of the magnetic stretch- 
ing against the growth of global magnetic helicity and magnetic 
helicity at the largest scales. In the time pictured, two shear- 
busts occur within 5 t^, and a clear double-peak structure is also 
visible in the magnetic helicity. 

The magnetic helicity grows particularly in low wavenum- 
bers k, associated with the growth of an isolated structure with a 
strong helicity polarity; this low-A: growth is a signature of an on- 
going inverse spectral transfer of magnetic helicity common for 
3D M HD systems ( Muller et al.|2012[ [Biska mp 2003; |Alexakis 
|et al.|2 008 ). The dramatic change in the bias of magnetic helicity 
in the system during one shear-burst is shown in Fig. l6| in Fig. 
Blthe large-scale magnetic helicity of the structures spawned has 
negative polarity in the two shear-bursts pictured. Large-scale 
magnetic helicity structures persist longer than the high-energy 
shear streams, and longer than it takes for the global energies to 
taper off. Because the magnetic helicity experiences an inverse 
cascade and our system has small dissipation this is theoretically 
expected. 



Shear-bursts change the magnetic helicity because they gen- 
erate significant currents through magnetic shear When no 
shear-burst is present in the system small filaments of high- 
current are common and likely indicate slow reconnection on 
small scales. However when a shear-burst grows, large-scale 
high-current structures grow at the same time. A typical growth 
in current around a shear-burst is shown in Fig. |7] 
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Fig. 6. (Above) probability distribution and (below) cumulative dis- 
tribution of local magnetic helicity, i.e. dH^ = A • B in the simulation 
volume before and during a shear-burst in simulation gl. During this 
shear-burst, as the global magnetic helicity grows, the tail of negative 
magnetic helicity grows, and grows much higher than the positive tail. 



1. Conclusions 

We have isolated a characteristic mechanism of dynamo ac- 
tion in MHD convection that operates through spontaneously- 
developing intermittent bursts of high shear. Because this pro- 
cess is robust and general, it is of potential importance for astro- 
physical small-scale dynamo action in turbulent convection sce- 
narios. These shear-bursts consist of the formation of large, co- 
herent, highly-sheared flows along-side magnetic structures with 
a magnetic -helicity polarity bias. The increasing shear causes 
a gradual rise of energy at all spatial scales due to magnetic 
stretching in the system over several buoyancy times. After some 
time, the shear flows become so highly aligned that their violent 
shearing breaks up the flow structure. Once the shear flows are 
destroyed, the elevated energy dissipates over several buoyancy 
times. Closely-spaced shear-bursts often occur in a series, cre- 
ating long periods of time where magnetic energy is elevated 
significantly above the steady-state. 
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Fig. 7. Current-squared in the lowest wavenumbers k = 1,2,3,4 
scaled to its initial time-value during a shear-burst in simulation g2. 
Global current- squared and magnetic shear are shown scaled and shifted 
to fit on the same scale for reference. 
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